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I. INTRODUCTION 



With the confirmation of atmospheric muon-neutrino oscillations by the Super-Kamiokande 
collaboration [lj], the need for physics beyond the standard model became clear. This prompted 
many authors to construct mass matrix models to explain the neutrino mixings and the still not 
well-determined neutrino mass spectrum. Meanwhile underground and reactor experiments P] 
confirmed the existence of solar electron-neutrino oscillations as a solution to the solar neutrino 
puzzle, whereby the detection of solar electron-type neutrinos was depleted [jj relative to the 
standard solar model projections [4]. Since then many of the first round of models have fallen 
by the wayside, as the once-preferred small mixing angle (SMA) solution for the solar neutrino 
oscillations has been replaced by the large mixing angle (LMA) solution. Accelerator experiments 
have now also contributed to the greatly increased precision of all observed oscillation results. 
In particular, the observed large atmospheric neutrino mixing has persisted as a nearly maximal 
Vu — v T mixing The solar neutrino mixing is now known to be large but not maximal [6]. On 
the other hand, only a relatively small upper limit has been placed on the v e mixing with the other 
two antineutrinos by two reactor neutrino experiments carried out several years ago [7j. Clearly 
the neutrino mixing pattern is totally unlike that observed in the quark sector. 

Despite the refinement of the experimental results, many of the more recent neutrino models 
have survived todate due to the uncertainty in the precise magnitude for the reactor angle, $13, 
and the unknown mass hierarchy for the neutrino spectrum: normal, degenerate, or inverted. In 
this paper we are primarily concerned with the model predictions for the reactor angle. While 
new reactor experiments Q are being planned or are already in construction to measure this angle 
down to 3°, i.e., sin 2 2#i3 ~ 0.01, it is of considerable interest to learn whether that reach will be 
sufficient to determine the angle through the detection of a v e disappearance signal, or whether 
a neutrino factory will be required which can probe a much smaller angle [jj]. For this purpose, 
we have surveyed 63 models in the literature which are still viable candidates and have reasonably 
well-defined predictions for #13 . Roughly half of the models predict that sin 2 covers the range 
from 0.015 to the present upper bound of 0.15 [7]. Hence half of the models can be eliminated 
in the next round of reactor experiments, based on the presence or absence of an observed v e 
disappearance signal. While none of the models proposed so far may be correct, the distribution 
of results for sin 2 2#i3 does provide some indication of what one may expect to find. 

In Sect. II we define the mixing angles and state the present experimental results. A brief 
description of the types of models based on lepton flavor symmetries and/or grand unification 
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for the quarks and leptons is presented in Sect. III. Several more comprehensive reviews of mass 
matrix models exist in the literature 1(| , and we refer the interested reader to them for additional 
details. We have restricted our attention to models with just three active neutrinos and no sterile 
neutrinos. This seems justified at present in light of the conflicting evidence for an additional 
oscillation between and v e observed by the LSND Collaboration but not by the Karmen 



Collaboration [2] . The MiniBooNE Collaboration is expected to throw more light on this situation 
shortly We make clear our acceptance criteria for the three flavor neutrino models to be 

considered and give tables of the mixing angles predicted for each in Sect. IV. Histograms for 
sin 2 #13 are plotted for all 63 models together, for the 27 models which give numerical predictions 
for all three mixing angles lying within the 2a experimental bounds, and separately for those 
models with normal or inverted neutrino mass hierarchy. Our conclusions are presented in Sect. V. 

II. LEPTONIC MIXING MATRIX AND PRESENT EXPERIMENTAL INFORMATION 

Here for completeness we present a brief description of the well-known neutrino mixing matrix 
and give an up-to-date summary of the constraints on the mixing angles. 

The light left-handed neutrino flavor states are related to the neutrino mass eigenstates by the 
linear combinations 

K) = ^2(U„ L ) ai \iyi), (1) 

i 

where a = e, fx, r and i = 1,2,3. If neutrinos are Majorana particles, the U UL transformation 
matrix is obtained by diagonalization of the effective left-handed Majorana neutrino mass matrix 
according to 

M diag = U T lMvUvl: (2) 

where M v is model-dependent and is typically constructed from some basic symmetry principle, or 



follows from the seesaw mechanism 1J] in grand unified (GUT) models. For models in which the 
light neutrinos are assumed to be Dirac particles, the U v l unitary transformation matrix can be 
obtained from the bi-unitary transformation which relates the flavor basis to the mass basis as in 

Mp* = Ut R M?U UL . (3) 

The light Dirac neutrino mass matrix, M®, is also constructed according to some symmetry prin- 
ciple. 
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In models based on some assumed lepton flavor symmetry, the charged lepton mass matrix is 
assumed to be diagonal in the lepton flavor basis with the charged lepton masses m e , and m T 
ordered along the diagonal; thus Ul l , the counterpart of U VL , is just the identity matrix. For GUT 
models, on the other hand, the charged lepton mass matrix is generally not diagonal in the GUT 
flavor basis. In this case in analogy with Eq. ([3]) above, Ul l can be determined from the bi-unitary 
transformation 

M d L m9 = UI R M L U LL , (4) 

where we have again adopted the convention that the right-handed fields act on the left and the 
left-handed fields on the right of the Dirac mass matrix. For either type of model, the left-handed 
neutrino Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is then given by [151 ] 

VpmNS = Ut U UL = Upmns®. (5) 



It is convenient to choose by convention 

UpMNS 



/ C12C13 S12C13 S 1S e lS \ 



-S12C23 - Ci 2 S23Sl3e t5 Ci 2 C 23 ~ Si 2 S 2 3S 13 e t5 S23C13 
\ S23S12 - C 12 C 2 3Sy3e l5 -C 12 S 23 - S 1 2C 2 3S 13 e lS C23C13 / 



(6) 



in analogy with the quark mixing matrix, along with the Majorana phase matrix, 

* = diag(e <xi , e iX2 , 1), (7) 

in terms of the three mixing angles, 9i 2 , #23 an d #13; the Dirac CP phase, 5; and the two Majorana 
phases, xi an d Xi- With Dirac neutrinos, one is free to make phase transformations on both Ul l 
and U Ul , so <E> is just the identity matrix. But in the Majorana neutrino case where U„ L is defined 
by Eq. (|2|), an arbitrary phase transformation is not possible when one demands real diagonal 
neutrino mass entries. Hence the presence of the Majorana phase matrix is required in order to 
adopt the convention for Upmns specified in Eq. (jfij). 

We now summarize the numerical information for these mixing parameters, as given by Maltoni, 
Schwetz, Tortola, and Valle 



la] in a recent updated global analysis which incorporates all the latest 



results cited in [3|-[3|. Within la accuracy, they found 





= (7.3- 


8.5) x 10~ 5 eV 2 


Am 2 31 


= (2.2- 


3.0) x 10~ 3 eV 2 


sin 2 6*12 


= 0.26- 


0.36, 


sin 2 9 23 


= 0.38- 


0.63, 


sin 2 6>i3 


< 0.025. 





(8) 



4 



No information exists for the Dirac or Majorana CP phases, for the neutrino mass hierarchy, or 
for the Dirac vs. Majorana nature of the neutrinos. Of special interest to us is the upper bound 
on sin 2 #13 which has mainly been determined by a non-observation of a depletion of the v e flux 
from the CHOOZ reactor [7]. In Sect. IV we shall turn to the model selection criterion we have 
used to extract the various model predictions. But first we give descriptions of the types of models 
that have been proposed. 



III. DESCRIPTIONS OF MODELS 



Most neutrino mass matrix models fall into two broad classes: those based on a lepton flavor 
symmetry which applies only to the charged lepton and light neutrino mass matrices and those 
based on some grand unification scheme which applies to both leptons and quarks. The latter 
class typically involves some family unification group which unites the quarks and leptons at the 
high gauge unification scale, while a flavor symmetry applying to the corresponding members of 
different families may or may not be invoked. Moreover, right-handed singlet neutrinos belong to 
irreducible representations and provide an exquisite way to give ultralight masses via the seesaw 



mechanism 



141 ] to the observed left-handed neutrinos. In the purely leptonic models some lepton 



flavor symmetry is generally considered, but the ultralight 0.001 — 0.1 eV neutrino mass scales 
typically remain a puzzle. Several symmetry schemes have been proposed within each of the two 
classes. Exceptional cases involve a model based on anarchy with no flavor symmetry and models 
with sequential right-handed neutrino dominance but no unification group identified. 

Of considerable interest is whether a given model exhibits a normal or inverted mass hierarchy. 
Models of either hierarchy can be found in several categories, while for others only one or the other 
is definitely preferred. We shall give broad and general descriptions here for each type of model 
considered. As cited earlier, we refer the reader to the comprehensive reviews [101 ] for more details. 



A. Anarchy a Model of Flavor with No Flavor Symmetry 

Contrary to our general study of mixing predictions of neutrino mass matrix models with some 
specified symmetry, we first consider a model based on neutrino mass anarchy, where the neutrino 
mass matrix is completely random as no flavor symmetry is specified [l^]. With this assumption 
of no fundamental distinction among the three flavors of neutrinos and hence no preferred basis 
for the neutrino states, one can simply examine statistically the neutrino mixing matrix. Initially 
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the authors of [17( first used a Monte Carlo analysis to test how many sets of mixing angles passed 
imposed cuts and then later applied the invariant Haar measure to the mixing angles. Still later 
de Gouvea and Murayama [3] refined the analysis by applying the Kolmogorov-Smirnov statistical 
test to the single sin 2 #13 variable to obtain a more precise lower bound on its expected value which 
is quoted in the first table. 



B. Models with Lepton Flavor Symmetries 

It is customary for models of this type to be formulated in the leptonic flavor basis for which the 
charged lepton mass matrix is diagonal. The differentiating feature of these models then resides 
solely in the light left-handed neutrino mass matrix. 



1. /i — r Interchange Symmetry and L e — L M — L T Conserved Flavor Symmetry 



The most general neutrino mass matrix exhibiting a — r symmetry is given by, 

' a b 6 \ 
M v = bed 
^b d c J 

In the restricted case where a, b < c ~ d, with M u having the texture 

/ a'e b'e b'e ^ 



0) 



\Z Am ' 2 atm 



b'e l + e -1 
y b'e -1 1 + e J 



(10) 



a normal mass hierarchy is obtained when a' and b' are of order one and the size of e is determined 



by the ratio, yAm^/Am^, multiplied by an order one coefficient which is a function of a' and 
b' . One finds Am 2 2 ~ m\, Am^ = m 2 - m\ > 0, sin 2 26 , 2 3 = 1, tan 2 6>i 2 ~ T^-fr, and sin6>i 3 = 0. 

On the other hand, for the simple form a = c = d = 0, the symmetry remains unbroken 
and an inverted mass hierarchy is obtained. In this case Am 2 2 = Am^ = m 2 , Am 2 2 = 0, and 
sin 2 2#23 = sin 2 2#i2 = 1 with sin 2 2#i3 = 0, as both the atmospheric and solar neutrino mixings 
are maximal, while the reactor neutrino mixing vanishes. Neither of these sets of predictions are 
observed experimentally, so soft symmetry breaking of either matrix texture must be introduced 
to obtain an acceptable model with the initial hierarchy unaltered. 
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The special inverted hierarchy texture case cited above actually exhibits an enhanced L e — L^ — 
L T symmetry which is more generally of the form |l9l | 

/o b V \ 



(11) 



6 
\V j 

Independent of the relative magnitudes of b and b' , this rank-2 matrix leads to a neutrino mass 
hierarchy which is inverted. Again soft symmetry breaking must be introduced in order for the 
model to be experimentally viable. Examples of these two lepton flavor symmetries are grouped 
together in the tables presented. 



2. S3 Lepton Flavor Symmetry 



In the case of S3 lepton flavor symmetry involving the permutation group of three flavors applied 
to both rows and columns of the neutrino mass matrix, the most general texture is 2(J 

' a b b^ 

M v = b a b , (12) 
\^b b a J 

in terms of two independent parameters, a and b. Alternatively, one can consider as a basis in 
the flavor space the unit mass matrix with l's down the diagonal and the democratic mass matrix 
with all unit elements. For the unit matrix, the mass spectrum is clearly degenerate and all mixing 
angles vanish. On the other hand, the rank-1 democratic matrix yields a normal mass hierarchy. 
Introduction of soft S^-breaking terms involving the parameter c as in 

( a b b ^ 
M v = b a-c b + c (13) 
yb b + c a — c j 

still respects the fx — r exchange symmetry and yields a neutrino mixing matrix which has the 
tribimaximal form suggested by Harrison, Perkins, and Scott [21 ] : 

1 2/^6 l/y/3 N 
Upmns= -1/V6 W3 l/y/2 , (14) 
\ - W6 1/V3 -1/^2 J 
corresponding to the mixing relations sin 2 2#23 = 1, sin 2 9 12 = 1/3, and sin 2 #13 = 0. These results 
are close to their experimental values, so only small corrections may be required. However, the 
diagonal charged lepton mass matrix clearly does not obey the S3 symmetry being considered. 



3. A± Lepton Flavor Symmetry 



The permutation group S4 of four objects was first considered by Ma and Rajasekaran 22j as a 
discrete flavor symmetry. Its non-Abelian subgroups are S3, D4 and A4. What makes the subgroup 
A4 of the twelve even permutations of S4 of particular interest is the fact that it is also the smallest 
discrete subgroup of SO(3) which has at least one three-dimensional representation. In fact, there 
are just four irreducible representations, one triplet and three singlets. The three lepton doublets 
can be placed in the 3 while the three right-handed charged leptons are each placed in one of 
the singlets of A±. With three Higgs doublets also transforming as a triplet, one can construct 
the charged lepton mass matrix. On the other hand, if the right-handed neutrinos are placed in a 
triplet representation, the Dirac and Majorana neutrino mass matrices can be generated with Higgs 
singlets. This application of A^ and other variations show that tribimaximal mixing of neutrinos 
can also be achieved while alleviating some of the problems of the S3 flavor symmetry. 

4- Other Lepton Flavor Symmetries 

Other lepton flavor symmetries have also been considered in the literature including, for exam- 
ple, SO(3) and SU(3). Another very popular starting point involves the arbitrary assignment of 



texture zeros in the light neutrino mass matrix 



231 ] . By doing so, one eliminates those assignments 



which do not yield a neutrino mixing matrix mimicing the nearly tribimaximal mixing form. The 
positions of the texture zeros may then point the way to some underlying flavor symmetry. Some 
of these models are included in our study. 

Various attempts have also been made to extend the lepton flavor symmetries proposed to the 
quark sector. But without a grand unification symmetry framework, they have met with rather 
mixed results. 

C. Sequential Right-Handed Neutrino Dominance 

Another class of models which do not neatly fit into the lepton flavor symmetry class or the grand 



241 ] . Here three right-handed 



unification class are models with right-handed neutrino dominance 
neutrinos are introduced which have a strong hierarchical mass spectrum. In the absence of any 
family gauge symmetry, one can still deduce that the light neutrino mass spectrum is controlled by 
the sequential dominance of the right-handed neutrinos, i.e., the mass of the heaviest left-handed 
neutrino is determined largely by the mass of the lightest right-handed neutrino, etc. With this 
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type of model, authors have shown that near tribimaximal mixings can also be obtained. 



D. Grand Unified Models 



An alternative approach is to start with a vertical family unification symmetry at some grand 
unification scale [25]. One can then try to impose a flavor symmetry on this structure which 
relates the corresponding members of each family; however, in many cases, one simply adopts 
an effective operator approach which arbitrarily assigns certain operators to each element of the 
mass matrices. Models formulated in this framework are obviously much more ambitious than the 
previous ones dealing only with the lepton sector. The mass and mixing results for the leptons 
are highly constrained by the input parameters introduced for the quark sector, since the Yukawa 
couplings apply to both quark and lepton Dirac mass matrices. Nevertheless, successful models 
exist in the literature which are still viable, even after tighter constraints on both the quark and 
lepton mixing parameters have been obtained experimentally. We briefly discuss models in the 
following categories. 



1. 577(5) and Flipped SU(5) Models 

Grand unification with a high scale SU(5) symmetry was first proposed by Georgi and Glashow 
26]. where the quarks, leptons and left-handed neutrinos can be placed into 10 and 5 representa- 
tions. With the appearance of neutrino oscillations it was then suggested to place the right-handed 
neutrinos into SU(5) singlets. The Higgs fields are conventionally placed in the 5 and 24 represen- 
tations. This group symmetry with a minimal Higgs structure proved to be less interesting when 
the limit on the proton decay lifetime increased several orders of magnitude above the predicted 
range of 10 29 — 10 30 years. 

An alternative procedure is to consider flipped SU(5) [3] in which the charged lepton and 
right-handed neutrino are interchanged in the 10 and 1 representations, as well as the conjugate 
up and conjugate down quarks in the 10 and 5 representations with respect to the usual SU(5) 
assignments. While models of this type have been pursued in the literature, none of them have 
numerical predictions for the reactor angle #13, so we do not elaborate on them further. 
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2. 50(10) Models with Higgs in High Rank Representations 



The SO(W) grand unification symmetry is an economical and attractive one 25[, for all sixteen 
left-handed quark and lepton fields and their left-handed conjugates fit neatly into one 16 repre- 
sentation per family. Many models exist in the literature which differ from one another by their 
Higgs representation assignments and flavor symmetry imposed, if any. To appreciate this, it is of 
interest to note the following decompositions of the direct product of representations: 

16 16 = 10 s 12O a 126 s , 

_ ( 15 ) 

16 16 = 10 45 210, 

where in the first product the 10 and 126 matrices are symmetric, while the 120 is antisymmetric. 



For the simplest Higgs structure [28| one can assume the presence of one 10//, one 126//, and 
one 126// representations, where the latter is needed to preserve a D-flat direction at the GUT 
scale in supersymmetric models. The 10// contains two Higgs doublets, which appear in the 5 + 5 



SU(5) decomposition, and contribute to the Dirac mass matrices, while the 126// contains a Higgs 
singlet which contributes to the right-handed Majorana neutrino mass matrix. Recently additional 
Higgs fields in the 120//, 210//, and 45// representations have also been considered by model 
builders in order to fine tune their predictions. The Dirac mass matrices will then have symmetric 
or antisymmetric textures, or even matrix elements with linear combinations of the two forms. The 
high ranks of these Higgs representations, rank-3, 4, and 5 are somewhat disfavored in the string 
theory framework 



With the vacuum expectation value of the 126// appearing near the GUT scale and giving 
massive entries to the right-handed Majorana mass matrix, the conventional type I seesaw mech- 
anism [14] provides a ready explanation for the ultralight left-handed neutrinos: 

M v = -M^M^M N , (16) 

where Mn is the Dirac neutrino mass matrix and M/j is the right-handed Majorana matrix in the 
basis convention that the conjugate left-handed fields appear on the left and the left-handed fields 
on the right of M^j. It is easy to show that models involving a type I seesaw yield a normal mass 
hierarchy for the light neutrinos 31]. 



Some authors also allow the possibility that Higgs triplet VEV's exist in the 126// and can 
give non-zero entries to the left-handed Majorana mass matrix, Mll- One then generates the light 
left-handed neutrinos through the type II (or mixed) seesaw mechanism [30] ]: 

M v = M LL - M^M^ l M N . (17) 
10 



In this case, the light neutrino mass hierarchy can be normal, degenerate, or inverted, the latter two 
occuring when the Mll contributions are comparable to or larger than the type I contributions. 

One possible breaking pattern of SO(10) down to the SM gauge group has the Pati-Salam 
group 



321 ]. Sf7(4) x SU(2)l x SU(2)r, as the intermediate gauge group, 

50(10) -> SU(4) x SU(2) L x SU(2) R -► SU(3) x SU(2) L x U(1) Y (18) 
This breaking can be achieved with a minimal Higgs content that has one 10 //, one 45 //, one 



54// and a conjugate pair of 126// © 126//. Due to the left-right symmetry, the resulting mass 
matrices are symmetric. In addition, if the minimal Higgs content described above is utilized, one 
has the following relations: the up type quark mass matrix and the Dirac neutrino mass matrix 
are identical, while the mass matrix of the down type quarks and that of the charged leptons are 
identical, up to some calculable Clebsch-Gordan coefficients which, when combined with the family 



symmetry, can be used to obtain the Georgi-Jarlskog relations 



33J required by phenomenology. 



The Majorana mass terms for the right-handed neutrinos arise from coupling to the 126//. These 
intra-family relations among the mass matrices greatly reduce the number of free parameters in 
the Yukawa sector, making these models very predictive. 



3. 5*0(10) Models with Lopsided Mass Matrices 

While the minimal Higgs models discussed above naturally preserve R-parity when the 126// 
develops a VEV and lepton number is violated by two units, they suffer from the disadvantage that 
they become non-perturbative above the GUT scale due to the high rank of the representations. 
Models which do not share this problem can be constructed by using lower rank Higgs represen- 
tations [34] including: 10//'s, one 45//, and one or two pairs of 16// - 16//. At the GUT scale, 
VEV's of the 45 and the SU(5) singlet parts of a 16// - 16// pair break the SO(W) symmetry to 
that of the standard model. Near GUT scale masses are generated for the right-handed neutrinos 
by pairs of 16// S'f7(5)-singlet VeV's which form an effective 126//. Due to the nature of these 
VEV's, lepton number is broken but only by one unit, so R-parity is broken. Hence it is neces- 
sary to introduce a matter parity in order to preserve the distinction between particles and their 
super-partners, unlike in the higher-rank Higgs type of models described above where R-parity is 
automatically preserved. 

Vacuum expectation values for the doublets in the 5(10//), 5(10//), and 5(16//) are then as- 
sumed to be generated at the electroweak scale. Since the surviving 5 Vd VEV is a linear combina- 
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tion of the VEVs of the 10// and 16//, one finds that tan [5 = v u /vd can be in the range of 5 - 55 
rather than simply 55 when the ratio involves just the two doublets from the 10 //. While the two 
doublets in the 5(10//) and 5(10//) Higgs representations contribute in a symmetric way to the 
ij components of the Dirac neutrino and up quark mass matrices, and to the charged lepton and 
down quark mass matrices, respectively, the doublet in the 5(16//) representation contributes only 
to the charged lepton and down quark mass matrices in a lopsided fashion. This follows because 
d,L and t° L lie in a 10(16), while d c L and £l he in a 5(16) matter representation. The complete 
Froggatt-Nielsen tree diagram then makes clear that if a charged lepton ij mass matrix element 
receives a large contribution while the transposed element ji vanishes, the opposite will be true for 
the down quark mass matrix. This lopsided behavior for the charged lepton mass matrix can lead 
to a large lepton flavor violation in ji — > e + 7, for example. The corresponding branching ratio for 
the higher rank Higgs models tends to be one or two orders of magnitude smaller. 

In order to obtain a successful SO (10) GUT model of either type, one must not only be able 
to generate appropriate neutrino masses and mixings, but the quark masses and CKM mixings 
for the quark sector must agree with the observed values after evolution downward from the GUT 
scale. This imposes considerably more constraints on the model than are present with the purely 
leptonic models discussed in part A. For either type of 50(10) model, the appearance of a PMNS 
neutrino mixing matrix close to the tribimaximal mixing form is usually regarded as accidental, 
rather than reflecting a symmetry inserted at the outset as in the purely leptonic flavor models. 

4- E e and E$ ® E$ Models 

Some authors have pursued models based on the exceptional E% gauge group. The matter fields 
of interest are placed in the 27 dimensional representation, while the Higgs fields are placed in a 
27//, 351//, and/or 351'// representations. One is then faced with the problem of making massive 
many of the extra fields which are present in such high dimensional representations. Some progress 
has been made, but none of the models have any firm numerical predictions for the neutrino mixing 
angles, so we do not consider them further. 

Even more ambitious models have attempted to deal with Es <8> Es grand unified models which 
naturally arise in the heterotic string theory. One of the E$s is assumed to break down to E& x 
SU(3), while the other represents a hidden symmetry. Many of these models, as well as the Eq 
models discussed in the previous paragraph, are formulated in five or six dimensions. Again no 
firm numerical predictions for the mixing angles have been obtained. 
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IV. RESULTS FOR THE MODELS SURVEYED 



In the previous Sect, we have presented broad general descriptions of models in the categories 
considered todate. Here we present results for the various models in the literature. We begin by 
defining our selection criteria for the three neutrino flavor models to be included in our survey. 
First of all we require that the models give the LMA solution for the solar neutrino oscillations and 
that firm and reasonably restrictive numerical predictions be given for the reactor sin 2 #13 mixing 
parameter. We do not require that the other two mixing angles or the mass squared differences 
be predicted, but all mixing angles for which information is given are listed in the tables. Models 
which are clearly in conflict with the present neutrino oscillation data are not considered, though 
we have not imposed an upper limit on the prediction of sin 2 #13. Many of the models have evolved 
with time and have been updated by their authors. As such we have generally listed only the latest 
published or archived version, except in cases where some important variation has provided two 
noticeably different results for the reactor angle. Only single references are given in alphabetical 
order to the accepted models in each category. The interested reader can readily track down earlier 
references to each model, if they exist, and can use them to learn the specific details of a certain 
model. Finally, we note that we have arbitrarily selected May 2006 as the cutoff date in accepting 
models for our compilations. 

In Table I we list 26 models cited in [18] and 3s| - [3 which exhibit one of the lepton flavor 
symmetry types: anarchy, L e — L^ — L T , S3 or S4, SO(3), and texture zeroes. Only six of them 
have firm predictions for all three mixing angles. The predictions for sin 2 #13 cover the full range 
of possibilities, from the present upper bound of 0.025 for the CHOOZ limit down to 10"° or less 
which would clearly require information from a neutrino factory to measure the result accurately. 
In Table II we list seven models cited in 



55] 



58j] based on sequential right-handed neutrino 



dominance in which no particular GUT model is assumed. All of these models have a normal 
hierarchy with a restricted range of predictions for the reactor angle, i.e., sin 2 #13 > 10 -3 . 

In Table III we list 24 S'O(IO) models cited in 5j| - [3] based on Higgs fields in the 10//, 



126//, 126//, and possibly 120// or 45// dimensional representations. As such, their mass matrix 
elements receive symmetric or antisymmetric contributions. Note that all models based on a type I 
seesaw mechanism have normal hierarchy, whereas those with inverted hierarchy would be highly 
unstable [31]. Only the type II seesaw models permit a stable inverted hierarchy depending upon 
the interplay of the type I seesaw and left-handed Majorana Mil contributions. Only three of the 
models predict a value of sin 2 6*13 ^ 10~ 3 . 
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Finally in Table IV are listed six 50(10) models 80J - 85f] involving Higgs fields in the 10#, 
16 h, 16/f, and 45 h dimensional representations which exhibit lopsided entries for the down quark 
and charged lepton mass matrices. Four of them predict values for sin 2 #13 < 3 x 10~ 3 . The other 
two have predictions so near to the CHOOZ bound, they are on the verge of being ruled out. 

In order to illustrate these results better visually, we have plotted histograms for the sin 2 #13 
predictions. In Fig. 1 we show the results for all 63 models. Only the results for the 27 models 
which predict all three mixing angles and which lie within the 2a bounds of Eq. (8) are plotted in 
Fig. 2. In Figs. 3a and b we separate the models into those with normal and inverted hierarchies, 
respectively. On the log(sin 2 #13) scale, we have divided each power of ten interval into three 
equal parts. For those models which give a range of values which occupy several intervals, we 
have rescaled them, so that each model has the same normalized area on each histogram for the 
predicted hierarchy. 

Clearly the majority of models prefer normal hierarchy, with the two main exceptions being 
models based on L e — — L T symmetry or those with well-designed texture zeros. Roughly half 
of the 63 models have sin 2 #13 > 0.004, or sin 2 29\3 > 0.015. This feature remains true even for the 
smaller number of 27 models which can accurately predict all three mixing angles. Since the next 
generation of reactor experiments is expected to reach values of order 0.01 for sin 2 2^13, we can 
expect that roughly half of the models will be eliminated based on the presence or absence of an 
observed v e disappearance signal. The possible need for a neutrino factory to reach even smaller 
values of the reactor neutrino mixing angle, if necessary, will then become apparent. On the other 
hand, if a disappearance signal is seen, and its value for sin 2 #13 can be well measured, the number 
of surviving models will be greatly reduced. Determination of the mass hierarchy will narrow the 
number down even further. 



V. CONCLUSION 



From our survey we found that the predictions for the angle #13 range from zero to the current 
experimental upper limit. For models based on GUT symmetries, normal mass hierarchy can 
be generated naturally. Inverted hierarchy may also be obtained in these models with a type-II 
seesaw, even though some fine-tuning is needed. Predictions for the mixing angle #13 in these 
models tend to be relatively large, with a median value sin 2 20\^ ~ 0.015. On the other hand, 
models based on leptonic symmetries can give rise to inverted mass hierarchy, and the predictions 
for #13 can be quite small. Therefore, if the inverted mass hierarchy is observed experimentally 
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and the mixing angle #13 turns out to be tiny, this experimental evidence will then give strong 
support to models based on lepton symmetries. However, if #13 turns out to be relatively large, 
one will not be able to tell the two different classes apart. A precise measurement for the deviation 
of 023 from 7r/4 can also be crucial for distinguishing different models. This is especially true for 
models based on lepton symmetries in which the deviation strongly depends on how the symmetry 
breaking is introduced into the models. Clearly precision measurements are indispensable in order 
to distinguish different classes of models and to narrow down the number of acceptable models. 
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TABLE I: Mixing Angles for Models with Lepton Flavor Symmetry. 



Reference 


Hierarchy 


sin 2 2# 23 


tan 2 9 12 


sin 2 13 


Anarchy M( 

dGM [18] 


>del: 

Either 






> 0.011 @ 2cr 


L e — — L r 
BM [35J 
BCM [36] 
GMN1 [37] 
GL [38] 
PR [39] 


Models: 

Inverted 
Inverted 
Inverted 
Inverted 
Inverted 




> 0.52 
< 0.58 


0.00029 
0.00063 
< 0.01 


> 0.007 


S 3 and S 4 rv 

CFM [40] 
HLM [41] 

KMM [42] 
MIN 1431 
MNY [44] 
MPR [45] 
RS [46] 

TY [47] 
T [48] 


lodels: 

Normal 
Normal 
Normal 
Inverted 
Normal 
Normal 
Normal 
Inverted 
Normal 
Inverted 
Normal 


1.0 
1.0 
1.0 

23 > 45° 
6 23 < 45° 
0.93 


0.43 
0.44 

0.43 


0.00006 - 0.001 
0.0044 
0.0034 
0.000012 
0.00z4 
0.000004 - 0.000036 
0.006 - 0.01 
< 0.02 


0.0025 
0.0016 - 0.0036 


A4 Tetrahec 

ABGMP [49] 
AKKL [50] 
Ma [51] 


iral Models: 

Normal 
Normal 
Normal 


0.997- 1.0 
1.0 


0.365 - 0.438 
0.45 


0.00069 - 0.0037 
0.006 - 0.04 



SO(3) Mode 

M [52] 


Is: 

Normal 


0.87- 1.0 


0.46 


0.00005 


Texture Zer 

CPP [53] 

WY [54] 


Models: 

Normal 
Inverted 
Inverted 
Either 
Either 
Either 






0.007 - 0.008 
> 0.00005 
> 0.032 
0.0006 - 0.003 
0.002 - 0.02 
0.02 - 0.15 
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TABLE II: Mixing Angles for Models with Sequential Right-Handed Neutrino Dominance. 



Reference 


Hierarchy 


sin 2 20 23 


tan 2 6 12 


sin 2 6»i3 


D [55] 


Normal 






0.008 - 0.14 


EH [56] 


Normal 


0.98 


0.32 


0.014 




Normal 


0.98 


0.34 


0.012 




Normal 


0.99 


0.45 


0.0009 




Normal 


0.97 


0.30 


0.014 


H [57] 


Normal 


1.0 


0.42 


0.0033 


K [58] 


Normal 


0.99 - 1.0 


0.40 - 0.62 


0.0027 



TABLE III: Mixing Angles for 50(10) Models with Symmetric/ Antisymmetric Contributions. 



Reference 


Hierarchy 


sin 2 2# 23 


tan 2 6 12 


sin 2 6>i3 


BaMa 


[59] 


Normal 


0.88 


0.33 


0.015 - 0.028 






Normal 


0.98 


0.44 


0.013 






Inverted 


0.88 


0.29 


0.024 


BMSV 


[60] 


Inverted 






> 0.01 


BKOT 


[61] 


Normal 


0.98 


0.28 


0.0001 - 0.0006 


BO 


[62] 


Normal 


0.98 - 1.0 


0.29 - 0.46 


0.0014 


BN 




Normal 


1.0 


0.36-0.39 


0.0009 - 0.016 


BeMa 


[64] 


Normal 


0.93 


0.40 


0.012 


BRT 


[65] 


Normal 


0.99 


0.35 


0.0024 


BW 


[66] 


Normal 






0(0.01) 


CM 


[67] 


Normal 


1.0 


0.41 


0.014 


DR 


[68] 


Normal 


0.98 


0.40 


0.0025 


DMM 


[69] 


Normal 






0.0036 - 0.012 


FO 


[2PJ 


Normal 


0.90 


0.31 


0.04 


GMN2 


[ZU 


Normal 


< 0.91 


> 0.52 


0.026 


KR 


[72] 


Normal 


0.93 


0.44 


0.058 


O 


[73] 


Normal 


0.94 


0.46 


0.0007 


Ra 


[74] 


Normal 






0(0.01) 


Ro 


[75] 


Normal 
Inverted 






0.0056 
0.036 


ST 


[76] 


Normal 


0.99 


0.46 


0.0001 - 0.04 


SP 


[77] 


Normal 


0.99 


0.42 


0.0002 


VR 


[78] 


Normal 


0.99 - 1.0 


0.40 - 0.61 


0.024 


YW 


[79] 


Normal 


0.96 


0.40 


0.04 
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TABLE IV: Mixing Angles for SO(10) Models with Lopsided Mass Matrices. 



Reference 


Hierarchy 


sin 2 2# 23 


tan 2 9 12 


sin 2 #i3 


A [80] 


Normal 


0.98 - 1.0 


0.38 - 0.50 


0.002 - 0.003 


AB [81] 


Normal 


0.99 


0.49 


0.0002 


BB [82] 


Normal 


0.97 


0.40 


0.0016 - 0.0025 


JLM [83] 


Normal 


1.0 


0.41 


0.019 


Mae [84] 


Normal 






0.048 


P [85] 


Normal 


0.99 


0.17 - 0.29 


0.0004 - 0.0025 
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Predictions of All 63 Models 
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FIG. 1: Histogram of the number of models for each sin #13 including all 63 models. 
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FIG. 2: Histogram of the number of models for each sin 2 #13 that give accurate predictions for all three 
leptonic mixing angles. 
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Models with Normal Hierarchy 
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FIG. 3: Histograms of the number of models for each sin #13 where the upper diagram includes models 
that predict normal mass hierarchy, while the lower diagram includes models that predict inverted mass 
hierarchy. 
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